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1 Introduction. 

In recently published papers [ 1,2 ] one of the authors has proposed a method to include the running 
coupling constant into the dynamics of rcggcizcd gluons based on the so-called "bootstrap condition" [ 
3 ]. This is a relation between the reggeized gluon trajectory 

(j{q) = -UsN^q^ / -TTT-^qlql, q = qi + q2 (1) 



and the gluon pair interaction given by the BFKL kernel [ 4 



K,{q„q[) = -2T,T,a, i (4 + 4)?^^ " ^ ) (2) 

where T is the gluon colour vector and qi + q2 = q[+ q2 = q- For the gluon channel T1T2 = —Nc/2 
and integrating the kernel one obtains the bootstrap relation [ 3 ] 

J {d'q[/{2nf)Kf-°^iq,q„q[)=L0{q)-ujiq,)-i0{q2) (3) 

This relation guarantees that in the t channel with the gluon colour quantum number the two gluon 
system has the reggeized gluon as its state. As a result the production amplitudes in the one-reggeized- 
gluon-exchange approximation, which serve as an input in the BFKL theory, satisiy unitarity in the 
leading order [ 5 ]. Thereby the whole scheme becomes self-consistent: otherwise one should add to the 
input amplitudes corrections following from the unitarity. Thus the bootstrap is a crucial element for 
the reggeization of the gluon and for the theory of the reggeized gluons as a whole. Therefore the only 
way to introduce the running coupling constant in a manner compatible with the gluon reggeization 
is to preserve the bootstrap. This was the motivation of the papers [ 1,2 ] 

Technically this can be achieved if one notices that (3) remains valid if both in the gluon trajectory 
and interaction each of the momenta squared is substituted by an arbitrary function of it. Indeed, if 
one takes for the trajectory 

and for the interaction kernel 

^qKQl, Ql) - --L1-L2 \ (—TTT + —r-rr)—. yr nTTTT (^j 

\ mi) viq2) viqi-Qi) vWi)v{q2)J 

then (3) will obviously be satisfied as well. The fixed coupling BFKL theory corresponds to a particular 
choice 

r^iq) = <?V(2as) (6) 

The idea of [1,2] was to change rj{q) so that it correspond to a running rather a fixed coupling. For 
the running coupling some conclusions about the form of rj{q) can be made considering the vacuum 
channel equation in the limiting case of very large q. Then putting (4) and (5) into the equation and 
assuming that r]{q) grows with q, in the leading approximation one gets an evolution equation 

dcf>{q\x)/dlnl/x={N,/27r)J ^^-/-C??,^) (7) 
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Comparison with the GLAP evolution equation in the leading order in In 1 jx (that is, in the double 
log approximation) allows to find the asymptotic form of 

~ gV(2a«(g2)), g ^ oo (8) 

It evidently differs from the fixed coupling case by changing the fixed coupling constant to a running 
one as(g'^). As a result, with a running coupling, both the gluon trajectory and its interaction have to 
be changed simultaneously in an interrelated manner, so that the resulting equation is different from 
the BFKL one already in the leading order. 

The behaviour of //(q') at small q, comparable or even smaller than the QCD parameter A, cannot 
be established from any theoretical calculation, since this domain is nonperturbative. In [ 1,2 ] these 
confinement effects were parametrizeed by an effective "gluon mass" m, choosing 77(g) in the form 

r?(g) = (&/2^)(<?' + m^) ln((g2 + m')lt^) (9) 

with h = (1/4)(11 — {2/3)Nf) and to > A, which agrees with (8) for large q and remains finite up to 
q = 0. 

A preliminary study of the properties of the pomeron with ri{q) given by (9) was performed by 
the variational technique in [ 1,2 ]. It was found that the intercept depended on the ratio m/A quite 
weakly: as m/A changes from 1.5 to 5.0 the intercept (minus one) A falls from 0.4 to 0.25. On the 
other hand, the slope depends on this ratio very strongly. This allows to fix the ratio to/ A to values 
in the interval 3.0-j-4.0. 

This variational study, although very simple, cannot however give values for the intercept and 
especially for the slope with some precision. Still less can be found by this method about the properties 
of the pomeron wave function essential for the high-energy behaviour of the physical amplitudes [ 6 
] . Finally, one does not receive any knowledge about the existence of other solutions with a positive 
intercept. All these reasons give us a motivation to undertake a numerical study of the two-gluon 
vacuum channel equation with the gluon tralectory and interaction given by (4) and (5) respectively 
and the function 77(g) satisying (8). Its concrete form has been chosen to be slightly more general 
than (9): 

77(g) = (6/27r)/(g) (10) 

where 

/(g) = (g2+TO?)ln((g2+m2)/A2) (11) 

It allows for the freezing of the coupling and the confinement proper to occur at somewhat different 
scales (to and toi respectively). However, on physical grounds, one feels that they should be of the 
same order. 

In sections 2 and 3 we present the basic equations in the form suitable for numerical analysis for 
the cases g = (forward scattering) and g 7^ 0. In Section 4 we describe the method of the solution 

and present the numerical results for the intercept, slope and the wave function at g = 0. The results 
for the intercept and slope, on the whole, agree with those found in [ 1,2 ] by the variational approach. 
An interesting new result is the existence of a second pomeron with the intercept, roughly speaking, 
two time less than for the leading one, but still positive. In Section 5 these results are applied to study 
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the asymptotical behaviour of the cross-section for the 7*7 scattering. Our conclusions are presented 
in Section 6. 



2 Basic equations. Pomeron at g = 

We consider the physical case Nc = 3. The units are chosen to have A = 1. In relating to observable 
quantities we take A = 0.2 GeV. 

The pomeron equation is the eigenvalue equation 

(-a;(9i) - coiq^mqi) + J [d^ /{2^fW,'^'{q^,(^^)<t>{q[) = E{q)(t>{qi) (12) 

where the "energy" eigenvalue E{q) is related to the pomeron trajectory via 

a{q) = 1 - E{q) = 1 + A - a'q^ (13) 

The last equation, valid for small, q defines the intercept A and the slope a' . In (12) the trajectories 
LO and the kernel JC""^ are given by the Eqs. (4) and (5) with T1T2 = —3 and the function r] given by 
(10) and (11). To symmetrize the kernel we pass to the function 

^(gi) = 4>{qi)Nv{qi)ri{q2) (14) 

We also take out the common numerical factor 6/ ((11 — 2/2>Np)'k) and express all terms via the 
function f{q) defined by (11). Then the equation for takes the form 

Md^mQ^) + j d^q'Mqu q[)i'{q[) = eiq^qi) (15) 
Here the "kinetic energy" is 

The interaction kernel consists of two parts, a quasilocal and a separable ones: 

L = + (17) 

They are given by 



and 



rW(, ,')- /fa) 1 //(92) //(g2) 1 f{q[). 

4'^''Hqi,q[)= , ■^^'^^ (19) 

x//(9i)/(<Z2)/(gi)/(g^) 

Both parts are evidently symmetric in qi and q[ . The scaled energy e is related to the initial one by 

Eq. (15) simplifies in the case when the total momentum of the two gluons is equal to zero. With 
q = the two parts of the kinetic energies become equal and the square roots in (18) turn to unity. 
So at g = the equation retains its form (15) with 

, . . f dXfjq,) 



5 



and the interaction given by (17) where now 
has really become local and 

This is the equation which we shall solve numerically. 

To make it one-dimensional we introduce the angular momentum of the gluons n and choose the 
solution in the form 

V-W = V'n(?^)expm0 (24) 

where ^ is the azimuthal angle. Integrating over it in the equation, we obtain an one-dimensional 
integral equation for the radial function ipniq^)' 

Ao{q)MQ^) + J dqlLn{q\ qlMqf) = e^„(g') (25) 
with the kernel now given by 

Ln{q^ ql) = -Bn{q\ql) + .5„o7r ./^^"^ . (26) 

where 

„ , 9 o, Z"^" ,, cosnd) 

Jo f{q^ +qt-2qqi cos (p) 

Note that Aq can be expressed via Bq: 

Ao{q) = (1/2) I dq!Bo{q',ql)j^^ (28) 

Evidently Eq. (25) is very similar to a two-dimensional Shroedinger equation with an attractive 
interaction provided by the local term and a positive kinetic energy described by A, which however 
grows very slowly at high momenta (as Inlng' according to [ 1 ]). Evidently the attraction becomes 
smaller with growing n. So wc expect to find negative energies, corresponding to intercepts larger 
than unity, only for small n. Remember that for the BFKL pomcron only the isotropic state with 
n — has a negative energy. Our calculations reveal that the introduction of the running coupling 
does not change this situation: states with \n\ > all have positive energies. So in the following we 
consider the case n = 0. 



3 Pomeron at g 7^ 0: the slope 

With q ^ the pomeron equation becomes essentially two dimensional. Rather than to attempt 
to solve it numerically at all q we limit ourselves to small values of q and determine not the whole 
trajectory a{q) but only the intercept a' defined by (13). This can be done in a much simpler manner 
using a perturbative approach. We present "the Hamiltonian" in (15) 

Hg=Ag + Lg {29) 

in the form 

Hg = Ho + W{q) (30) 
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and calculate analytically W{q) up to terms of the second order in q. Then for small q the value of 
the energy e{q) will be given by the standard perturbation formula 

e{q) = e(0)+ < W{q) > (31) 

where < > means taking the average with the wave function at g = 0, determined from the numerical 
solution of the equation discussed in the previous section. Thus we evade solving the two-dimensional 
problem, but, of course, cannot determine more than the slope. Fortunately it is practically all we 
need to study the high-energy asymptotics (although, of course, the knowledge of the trajectory as a 
whole might be of some interest). 

In order to derive an expression for W{q) we pass to the relative momenta / and I' 

gi = (l/2)g + l, q2 = {l/2)q-l, q[ = {l/2)q + I' , q'^ = {l/2)q - I' (32) 

Up to the second order in q we have 

f{q,) = /(0[1 + ai(g ■ + '^q' + y (g • if] (33) 

where 

a, = a,il) = [1 + ln(/^ + m^) - ^^^] (34) 
r 1 — m? T 1 

«^ = «^ W = [pT^ + jpi^] 7(^y (35) 

The expansion for f{q2) differs by changing the sign of I (or of q); for f{q[) and f{q2) it suffices 
to replace I by I' in the expressions for f{qi) and f{q2)- We use also the notation a'l = ai(Z') and 
O2 = a2{l'). We also need the expansion for f{q): 

f{q)=fmi + asq') + 0{q') (36) 

where 

The perturbation W{q), up to second order in q, can be expressed via the introduced functions ai,2 
and a'12 and the constant a^. After some calculations we find a part of W coming from the kinetic 
term in Hamiltonian in the form 

The part of W coming from the quasilocal part of the interaction can be written as 

W2{l,l') = l[ai{q-l) - a[{q-n]^Li''\l,n (39) 
and the one coming from the separable part as 

Ws{l,l') = [{as - ^l±^)q' - lia2 - al){q ■ if - ^(4 - a'l){q ■ l'f]Lt'\l,l') (40) 
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As mentioned, only isotropic solutions have the intercept larger than one and are of interest. 
Then the expression for W{q) = Xli=i 2 3 ^» *° integrated over the azimuthal angles. Thus 
intergrated values will be denoted Wj, i = 1, 2, 3. Using (26)-(28), they can be conveniently expressed 
via the kernel -B„ (eq. (27)): 

i^ = \J ^^";^{[-T + - y)«'>o(/,/') - a,a[ll'B,il,n} + Ml^a,+a,P) (41) 
^ = + all'^)Bo{l, I') + a,a[ll'B^{l, I') (42) 

The slope is given by the momentum average of the sum of these expressions, taken with a given 
isotropic wave function: 

, ^ (1/2) / diM"mm{w, +ws)+j dmrw, 

4 Numerical procedure and results 

Eq. (25) was first changed to the variable t = Ing^ whereupon the wave function and the kernel 
transform according to 

V'(g') ^ m = #(g2) (45) 

and 

L{q\ql) L{t,t^) = qq,L{q^,ql) (46) 

Then the equation was reduced to a finite system of linear equations by approximating the integral 
by a sum 

dtF{t)o^Y.WiF{ti) (47) 

with points U and weights Wi depending on the chosen approximation scheme. The final equation is 
thus 

n 

BijXj = exi, i = I, ...n (48) 

where 

x^ = V^V'(ii) (49) 

and 

Bij = A{ti)6ij + ^WiWjL{ti, tj) (50) 

A somewhat delicate part of procedure proved to be the integration over the angle in (26), since at 
large values of q^ and qf the integrand becomes strongly peaked at = 0, so that one should take 
much care to obtain a reasonable precision. 

After determining the lowest eigenvalues e and the corresponding eigenvectors Xi the wave function 
in the momentum space is directly given by (45) and (49) at points q^ = exp ij. It should be normalized 
according to 



f 

J — ( 



d^q 
(27r)2 



\i>{i)V = 1 (51) 
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Note that this wave function is a partially amputated one (see Eq. (14)). The full (nonamputated) 
wave function is given by $((7) = ip{q) /r]{q). It is this function that appears in the physical amplitudes. 

The results for the lowest (and negative) eigenvalues of energy for the case n = (isotropic 
pomeron) are presented in Figs. 1 and 2. Actually energies with an opposite sign are shown, which 
according to (13) give precisely the intercepts (minus one). As mentioned, the QCD scale here and 
in the following is taken to be A = 0.2 GeV. In Fig.l the intercepts are shown for the case when the 
two scales m and mi in (11) are equal. Fig. 2 illustrates the dependence of the intercepts on the 
ratio m/mi. The most interesting observation which follows from these figures at once is that in all 
cases one observes two positive intercepts, which correspond to two different supercritical pomerons, 
the leading and subleading ones. The intercept of the leading pomeron is found to be in accordance 
with our earlier calculations, performed by the variational method (it is slightly larger, which was to 
be expected). For physically realistic values of m and nii in the interval 0.5 1.0 GeV it takes on 
values in the region of 0.5 0.3 falling with the masses m and mi. The same trend is seen for the 
intercepts of the subleading pomeron, which lie in the interval 0.25 0.15. 

The slopes of the two found pomerons are shown in Fig. 3 as a function of m for the case m = mi. 
The slopes depend very strongly on the value of the regulator mass. The physically reasonable slopes 
for the domimant pomeron of the order of a' 0.25 {GeV/c)~'^ restrict the values of m to the region 
0.7 - 0.9 GeV. So finally we choose 

m = 0.82 GeV (52) 
which leads to the following parameters of the leading (0) and subleading (1) pomerons 

Ao = 0.384, a'o = 0.250 (GeV/c)-'^; Ai = 0.192, a[ = 0.124 (GeV/c)-'^ (53) 

In Fig. 4 we show the coordinate space wave functions $(r) of these two pomerons. 



5 Pomerons and the high-energy scattering. Discussion 

To apply the found results to the actual physical processes one has to couple the pomerons to the 
external sources corresponding to the colliding particles. The only way to do it in a more or less 
rigorous manner is to assume that both the projectile and target are highly virtual photons with 
momenta q and p respectively, —q^ — >> A^ and —p^ ~ >> A^. Then the nonperturbative 
effects inside the target and projectile can safely be neglected. As shown in [6] the scattering amplitude 
in the high colour number limit then takes an eikonal form for fixed transverse dimensions of the 
projectile and target and leads to a cross-section 

cr = 2 y (fR(fr(fr'pg{r)pp{r') (1 - exp(-^(i/, R, r, r'))) (54) 

where 

/cPq(PqicPq' 
(27r)6 ^ ^' e^P*9^ n ~ expiqir){l - expiq^r') (55) 

is essentially the Fourier transform of the (nonamputated) Green function of Eq. (12), G{i^, q, qi, q[), 
considered as a function of the energetic variable u = pq and with q = qi+q2 = 9i + ^2 • The functions 
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Pq and pp correspond to the colour densities of the projectile and target photons, respectively. Their 
explicit form was found in [7] for both transverse and longitudinal photons. 

The found supercritical pomerons represent a part of the total pomeron spectrum which contributes 
to the Green function a term dominating at high energies 

Gp{iy,q,qi,q[) = ^ g2)$:(gi, 9^) (56) 

i=0,l 

where Ui and $i are the trajectories and wave functions of the leading (0) and subleading (1) pomerons. 
At high f we can neglect the dependence on the total momentum q of the wave functions, taking them 
at g = 0, and approximate the trajectories according to (13). Then all the quantities in (54) become 
determined, so that we can calculate the cross-sections for both the transversal and longitudinal 
projectile photon and thus find the structure function of the virtual photon target. We have taken 
for the latter a transversal photon with the lowest momentum admissible of P = 1 GeV/c. The 
resulting structure functions are presented in Fig. 5 for the interval of small x which we extended to 
extraordinary small values to clearly see the unitarization effects. 

To move closer to reality one has to consider hadronic target and projectiles. The confinement 
effects make any rigorous treatment of such a case impossible. Rather than to introduce arbitrary 
parameters (in fact, functions) we extend the formula (54) to hadronic target and projectile subsi- 
tuting the photonic colour densities by hadronic ones. For the latter we take a Gaussian form and 
a normalization which follows from the comparison to the electromagnetic densities with only the 
simplest quark diagrams taken into account. In particular for the proton wc take the Gaussian p, 
with the observed electromagnetic proton radius and normalized to three. Such a treatment, in all 
probability, somewhat underestimates the density, since it does not include coupling to gluons. 

The proton structure functions and the proton-proton total cross-sections which follow from this 
approximation for the densities are shown in Figs. 6 and 7 respectively. To see the unitarization 
effects we had again to consider extraordinary high values of 1/a; and energies, well beyond our 
present experimental possibilities. 

In discussing these results, we have first to note that their overall normalization is somewhat 
undetermined, since the exact scale at which u enters into In u factors is unknown. A second point to 
note is that the subleading pomeron contribution is always very small: it amounts to a few percent 
at smallest values of 1/x and s considered and naturally gets still smaller at higher 1/x or s. 

As one observes from Figs. 5-7, the structure functions and cross-sections monotonously rise with 
1/x, s and . Studying the asymptotics of the solutions of Eq. (12) at high q and of Eq. (54) one 
can show that this rise is logarithmic. In particular, the structure function of the virtual photon rises 
as ln^(l/x) and as In'^(g^) with /3 ^ 2.5. The proton-proton cross-section eventually rise as In^ s, 
as expected. Comparison to the Froissart bound (dash-dotted line in Fig. 7) shows however that it 
remains far from being saturated. 

The most interesting result that follows from Figs. 5-7 is that the unitarization effects become 
visible only at exceedingly very small values of x or very large values of s, well outside the range of 
the present experiment. They appear earlier at lower Q^. Still at the smallest value Q = 2 GeV/c 
considered, the exchange of more than one pomeron achieves only 15% of the total for the proton 
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structure function at a; = 10~^°. Likewise the relative contribution of many pomerons to the proton- 
proton cross-section rises to 23% only at s lO^Ge^. 

Comparing the calculated proton structure functions and the cross-sections with the experimen- 
tal results at highest 1/x and s available we observe that our results are essentially smaller than 
the observed ones. Experimental value of F'2p{Q'^,x) at Q"^ = 8.5 (GeVjcf and x = 0.000178 is 
1.19±0.05±0.16 [8 ]. Our calculations only give a value 0.17. The pp cross-section at ^/s = 1800 GeV 
is around 80 mbn [9], whereas our result is 18.5 mbn. Of course, having in mind the uncertainties in the 
overall normalization and a very crude picture for the pomeron coupling to the proton assumed, one 
cannot ascribe too much importance to this fact. However one is tempted to explain this underesti- 
mation of the experimental values by the simple reason that we are too far from the pure asymptotical 
regime yet and that other solutions of Eq. (13) different from the found supercritical pomerons and 
having their intercepts around unity give the bulk of the contribution at present energies. This may 
also explain the notorious discrepancy between a high value of the hard pomeron intercept, of the 
order 0.35-0.5, and the observed slow growth of the experimental cross-section, well described by the 
"soft pomeron" with an intercept around 0.08. 

If this picture is correct then we may expect that with the growth of energy the cross-sections will 
grow faster and faster, until at ^/s 10 TeV they will become well described by a pure hard pomeron 
with the intercept 0.35-0.5. This power growth will continue until energies of an order 1000 TeV 
when finally the unitarity corrections set in to moderate the growth in accordance with the Froissart 
bound. 

6 Conclusions 

The result of our study show that with the running coupling included the pomeron equation posseses 
bound state solutions which have negative energy and thus intercepts greater than unity. These 
solutions correspond to supercritical pomerons in the old sense, that is, they represent simple poles 
in the complex angular momentum plane. A new result is that two such pomerons appear. However 
the subdominant pomeron does not seem to play any role in describing the asymptotical behaviour of 
the amplitudes. The intercepts of the pomerons depend weakly on the infrared regulator parameter 
and stay in the region 0.35-0.5 for its physically reasonable values. The introduction of the running 
coupling and thus a scale provides for a notrivial slope for the pomeron, which is responsible for the 
physically reasonable behaviour of the cross-sections at very high energies. 

For realistic photonic cross-sections and with a rather crude approximation for the hadronic ones 
unitarization effects begin to be felt at extraordinary high energies, of the order 100 — 1000 TeV (or 
equivalently x < 10~^° — 10~^^). Until these energies a single pomeron exchange remains a very good 
approximation to the asymptotic amplitude. 

Comparison to the experimental cross-sections and structure functions at the highest energy (low- 
est x) achieved seems to confirm the widespread opinion that we are still rather far from the asymp- 
totical regime and that other states, different from the supercritical pomerons, give the dominant 
contribution. 
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9 Figure captions 

Fig. 1 Pomeron intercepts as a function of the infrared regulator mass m = mi; A = 0.2 GeV. 

Fig. 2. Pomeron intercepts for different values of the confinement parameter mi and the couphng 

freezing parameter m. 

Fig. 3. Pomeron slopes as a fimction of the infrared regulator mass m = mi; A = 0.2 GeV. 

Fig. 4. Coordinate space wave functions for the leading (^ol^")) cind subleading ($i(r)) pomerons. 

Both r and $ arc in units 1 /A ~ 1 fm 

Fig. 5. Structure functions of a virtual photon (P = 1 GeV/c) at = 4 and 100 {GeV/c)^ as a 
function of x (solid curves). Dashed curves show the contribution of a single pomeron exchange. 
Fig. 6. Proton structure functions at = 4 and 100 (GeV/c)'^ as a function of x (solid curves). 
Dashed curves show the contribution of a single pomeron exchange. 

Fig. 7. Proton-proton total cross-sections as a function of cm. energy (the solid curve). The 
dashed curve shows the contribution of a single pomeron exchange. The dash-dotted curve marks the 
Proissart bound. 
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